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商 要 : 冰川 物质 平衡 是 反映 气候 变化 的 重要 指标 ,对 于 区 域 生 态 环境 评估 和 灾害 防治 具有 重要 意 


义 。 采 用 Landsat 系列 遥感 
区 内 的 冰川 边界 ， p 


InSAR) 技 术 监 测 区 域 冰 川 形变 特征 ,并 分 析 冰 川 物质 平衡 过 程 。 
续 退 缩 , 且 近 10a 退 缩 趋势 更 为 显著 。 保 护 区 内 冰川 总 
变化 率 为 - el i 
围 内 冰川 退缩 最 为 显著 。(3) 2020 年 珠 峰 自然 保护 区 冰川 形变 速率 介 
形变 在 海拔 4200~4400 四 和 40"~45" 处 沉降 最 为 严重 。 
(4) 气温 上 升 、 降 水 减少 可 能 是 导致 珠 峰 自然 保护 区 冰川 物质 亏损 的 主要 因素 。 同 时 , 空 


峰 自然 保护 区 冰川 持 


5400~5600 m 和 10°~15° 范 
于 -129.069~140.252 mm*a 之 间 


SKIN 


感 影像 ,运用 比值 闵 值 法 和 和 目 视 解 译 法 ,获取 1990 一 2020 年 珠 峰 自 然 保护 
积分 布 与 变化 特征 ,同时 基于 差分 干涉 测量 短 基线 时 序 分 析 (SBAS- 


a ae mona 
只 退缩 247.16 km?, 
围 内 ,其 中 


x le] RE 


异 和 地 形 等 因素 也 可 能 是 导致 冰川 物质 平衡 差异 的 重要 因素 。 
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冰川 变化 与 气候 变化 息息相关 ,对 气候 系统 有 


重要 影响 "。 人 研究 表明 冰川 变化 对 温度 响应 程度 较 
高 ,其 变化 特征 能 反映 出 全 球 气候 变化 特点 ”。 随 


着 全 球 气候 变 暖 , 尤 其 是 21 世纪 以 来 ,冰川 退缩 趋 
势 愈 加 显著 。 一 定时 期 内 冰川 消融 有 利于 水 资源 
补给 ,而 长 期 的 冰川 消融 不 仅 会 引发 冰 湖 演 决 vk 
川 泥石流 等 自然 灾害 ,还 会 导致 海平 面 上 升 ,威胁 
沿海 居民 的 生命 安全 和 社会 经 济 发 展 。 因 此 ,研究 
全 球 变化 背景 下 的 冰川 变化 ,对 认识 区 域 冰川 变化 
及 其 影响 .评价 区 域 的 生态 环境 .探究 冰川 与 气候 
的 相互 作用 关系 具有 重要 科学 与 社会 意义 ,冰川 面 
积 的 分 布 与 变化 特征 也 一 直 是 科学 研究 热点 “1。 

尽管 冰川 面积 分 布 与 变化 特征 研究 对 冰川 区 的 灾 
害 预测 和 防治 .水 资源 利用 等 具有 重要 意义 ,但 物 
质 平 衡 作 为 冰川 变化 重要 指标 ,可 以 直观 反映 区 域 
气候 变化 特征 ,也 是 连接 冰川 与 水 资源 的 重要 枢 
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SAL! 。 因 此 ,研究 冰川 面积 的 同时 监测 冰川 物质 平 
衡 , 有 助 于 深化 对 冰川 响应 气候 变化 的 科学 认识 ， 
增加 对 冰冻 圈 生 态 系 统 演变 的 理解 ,为 全 球 气候 变 
暧 下 的 环境 应 对 策略 提供 科学 参考 "1。 

传统 的 冰川 监测 以 实地 考察 为 主 ,精度 高 ,但 
耗 时 耗 力 ,难以 实现 大 范围 动态 监测 。 和 遥感 技术 的 
发 展 使 冰川 研究 进入 了 朵 新 阶段 ,光学 遥感 数据 因 
其 直观 清晰、 易于 判读 的 优点 被 广泛 应 用 于 冰川 
面积 研究 中 。 而 现 阶段 冰川 表面 物质 平衡 的 变化 
监测 ,大 多 采用 物质 平衡 模型 和 多 源 数字 高 程 模型 
(DEM) 空 间 匹 配 的 大 地 测量 法 。 尽 管 物质 平衡 模 
型 在 冰川 物质 平衡 研究 中 发 挥 了 重要 作用 ,但 模型 
存在 参数 不 确定 性 和 区 域 适用 性 问题 。 目 前 ,多 源 
DEM 空间 匹配 的 大 地 测量 法 发 展 成 熟 且 被 广泛 应 
用 ,但 存在 多 源 数据 匹配 问题 。 冰 川 表面 形变 作为 
冰川 物质 平衡 的 重要 参数 ,可 直接 参与 反 演 冰 川 物 


基金 项 目 : 国家 自然 科学 基金 项 目 (42071277); 重 庆 市 教委 科学 技术 研究 计划 项 目 (KJQN202200526); 重 庆 市 教育 委员 会 人 文 社会 科学 
人 研究 项 目 (22SKSZ030); 重 庆 师 范 大 学 基金 项 目 (22XLB003); 国 家 级 大 学 生 创新 创业 训练 计划 项 目 (202210637008) 资 助 


作者 简介 : 费 琴 (1987-), 女 ,博士 ,副教授 ,主要 从 事 3S 技 术 应 用 


通讯 作者 : ; 


的 研究. E-mail: yunngy@126.com 
梁 文 莉 (1995-), 女 ,人 硕士 ,主要 从 事 鹏 感应 用 研究 . E-mail: liangwenli_ly@163.com 


1592 FA wv 46 卷 


质 平 衡 计 算 ,其 变化 趋势 能 直观 反映 冰川 物质 平衡 
变化 。 因 此 ,研究 冰川 表面 形变 特征 可 在 一 定 程度 
反映 冰川 物质 平衡 变化 。 近 年 来 ,雷达 遥感 尤其 是 
差分 干涉 测量 短 基线 时 序 分 析 (SBAS-InSAR ) 技 术 
已 经 在 毫米 级 地 表 形 变 应 用 中 取得 了 多 项 标志 性 
成 果 * ,其 高 精度 观测 、 较 强 适用 范围 以 及 可 长 时 
间 监 测 缓 慢 地 表 形 变 的 优势 中 ,为 冰川 表面 物质 平 
衔 研 究 提供 了 新 的 途径 。 

本 文 以 喜马拉雅 山 珠 峰 自然 保护 区 冰川 为 研 
究 对 象 , 采 用 Landsat 系列 光学 影像 和 哨兵 1 号 (Sen- 
tinel-1) 雷 达 数 据 , 运 用 比值 阔 值 法 和 SBAS-InSAR 
技术 ,监测 珠 峰 地 区 冰川 面积 分 布 .变化 特征 以 及 
冰川 表面 形变 ,分 析 冰 川 表面 形变 的 时 空 分 异 特 
征 ,同时 结合 前 人 研究 成 果 ,探究 冰川 物质 平衡 与 
气候 变化 的 响应 关系 ,以 期 为 珠 峰 保护 区 冰川 灾害 
的 防治 与 治理 提供 科学 依据 。 


1 研究 区 概况 


珠 峰 自然 保护 区 地 处 中 国 和 尼泊尔 交界 处 (图 
1) ,总 面积 约 3.38x10: km ,是 世界 上 海拔 最 高 的 自 
然 保护 区 """。 据 中 国 第 二 次 冰川 编目 显示 中, 保护 
区 内 共有 冰川 1049 条 ,总 面积 约 1450.40 km ,其 中 
绒布 冰川 分 布 面积 最 大 。 人 研究 区 为 大 陆 性 高 原 气 
候 , 光 照 充 足 , 冬 寒 夏 凉 。 嘎 尔 气 象 站 资料 显示 ,区 


域 年 均 气 温 2.1 % ,日 照 百分率 75.3% ,年 均 降 雨量 
270.5 mm, 两 热 同 季 。 珠 峰 自 然 保 护 区 因 其 独特 复 
林 的 地 理 位 置 和 气候 条 件 ,成 为 科学 研究 关注 的 热 
点 区 域 之 一 。 


2 数据 与 方法 


2.1 数据 来 源 
211 ZRP 本 文 使 用 的 Landsat 系 列 遥 感 影像 
来 源 于 美国 地 质 调查 局 (USGS) (https://earthexplor- 
er.usgs.gov/) 和 地 理 空间 数据 云 (https://www.gscloud. 
cn/)。 为 提高 冰川 边界 提取 的 精确 度 , 尽 量 选取 少 
云 . 少 积 雪 的 遥感 影像 ,最 终 使 用 的 影像 如 表 1 所 示 。 
Sentinel-1 卫星 是 欧洲 航天 局 为 监测 全 球 尺 度 
陆地 海洋 地 形 和 极地 等 设计 的 双星 观测 卫星 , 搭 
载 能 提供 连续 图 像 的 C 波 段 合成 孔径 雷达 。 其 升 轨 
影像 集中 在 夜间 成 像 , 降 轨 影像 在 白天 成 像 。 由 于 
冰川 热 融 可 能 导致 影像 失 相关 , 故 选择 更 稳定 的 升 
轨 数 据 ""。 本 研究 选取 2020 年 1 月 1 日 到 12 月 26 
日 共 72 景 升 轨 数据 ( 表 2) ,轨道 文件 采用 欧 空 局 提 
供 的 精确 轨道 数据 (https:We4ft01.cr.usgs.gov/) 以 消 
除 大 气相 位 影响 和 平地 效应 ,轨道 数据 精度 在 
10 cm 以内。 
2.1.2 DEM DEM 数据 采用 ASTER GDEM 和 
SRTM DEM 2 种 ,分 辨 率 均 为 30 m。2 种 DEM 均 可 


SN | Mo 


aa 


图 1 研究 区 概况 
Fig. 1 Overview of the study area 
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R1 遥感 影像 数据 


Tab. 1 Remote sensing image data 


年 份 行 号 / 列 号 获取 日 期 (年 -月 -日 ) 云 量 /% 卫星 /传感器 

1990 142/40 1991-10-12 1.00 Landsat4—5/TM 
141/40 1988-12-15 1.00 Landsat4—5/TM 
140/40 1989-11-09 1.00 Landsat4—5/TM 
140/41 1989-11-09 5.00 Landsat4—5/TM 
139/40 1990-06-14 1.00 Landsat4—5/TM 
139/41 1990-06-14 33.00 Landsat4—5/TM 

2000 142/40 2000-12-15 1.00 Landsat7/ETM+ 
141/40 2000-11-22 1.00 Landsat7/ETM+ 
140/40 2000-11-15 1.00 Landsat7/ETM+ 
140/41 1999-04-27 17.00 Landsat4—5/TM 
139/40 2001-12-29 0.00 Landsat7/ETM+ 
139/41 2000-12-26 2.00 Landsat7/ETM+ 

2010 142/40 2009-09-27 14.00 Landsat4—5/TM 
141/40 2010-12-12 1.00 Landsat4—5/TM 
140/40 2009-11-08 1.00 Landsat7/ETM+ 
140/41 2012-12-02 4.00 Landsat7/ETM+ 
139/40 2010-04-18 4.00 Landsat4—5/TM 
139/41 2012-10-08 5.00 Landsat7/ETM+ 

2020 142/40 2020-10-11 2.43 Landsat8/OLI_TIRS 
141/40 2020-10-20 4.24 Landsat8/OLI_TIRS 
140/40 2020-10-29 1.22 Landsat8/OLI_TIRS 
140/41 2020-10-29 3.65 Landsat8/OLI_TIRS 
139/40 2020-10-22 0.21 Landsat8/OLI_TIRS 
139/41 2020-10-22 10.62 Landsat8/OLI_TIRS 

表 2 Sentinel-1A 数据 SAR 技术 参考 DEM。 


Tab. 2 Sentinel-1A data 


参数 Sentinel-1 
波段 C 
产品 类 型 SLC 
成 像 模 式 IW Mode 
极 化 方式 VV 
轨道 方向 升 轨 
距离 向 分 辩 率 /m 5 
方位 向 分 辨 率 /m 20 


影像 开始 到 结束 日 期 2020 年 1 月 1 日 一 2020 年 12 月 26 日 


用 于 提取 冰川 边界 ,但 SRTM 受 镜面 反射 . 回 波 滞后 
和 雷达 阴影 影响 ,数据 会 出 现 缺 失 和 异常 ”。 而 
ASTER GDEM 数据 缺失 较 小 ,可 较 好 反映 地 貌 特 
征 。 此 外 , 龙 四 春 等 5 基于 不 同 参 考 DEM 进行 D- 
InSAR 技术 处 理 , 结 果 表 明 SRTM 的 差分 干涉 图 的 
噪声 较 少 。 故 本 研究 将 ASTER GDEM 用 于 冰川 边 
界 识 别 以 及 海拔 和 坡度 计算 ,将 SRTM DEM {EX In- 


2.2 研究 方法 
2.2.1 冰川 边界 提取 本 研究 选用 较为 精准 有 效 的 
比值 阐 值 法 (band3/band5) ,提取 珠 峰 自然 保护 区 的 
冰川 边界 。 对 比 多 次 实验 结果 ,确定 TMETM+ 影 像 
比值 图 像 的 最 佳 阔 值 为 1.8,OLI 影 像 为 1.0。 基 于 决 
策 树 分 类 得 到 二 值 图 像 ,最 后 进行 目 视 解 译 。 同 
时 ,基于 第 二 次 冰川 编目 数据 对 结果 进行 修正 ,从 
而 提高 冰川 边界 解 译 精 度 ( 图 2)。 
2.2.2 冰川 总 面积 计算 采用 年 均 变 化 率 评 估 珠 峰 
自然 保护 区 冰川 总 面积 变化 情况 ,计算 公式 如 下 : 

(AS,/S) 

T, 


APAC,= a= 100% (1) 


式 中 :APAC, 为 年 均 变化 率 ; AS, 为 第 i 时 段 内 冰川 
面积 变化 量 ; S 为 第 i 时 段 内 初始 年 份 冰 川 面积 ; 
AT, 为 第 ;时 段 的 时 间 间 隔 。 

2.2.3 冰川 形变 速率 监测 Berardino 等 "于 2002 年 
提出 小 基线 集 技术 (SBAS-InSAR 技术 ) ,该 方法 无 需 
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考虑 主 影像 选择 等 问题 ,对 时 空 失 相 干 和 大 气 等 
因素 不 敏感 , 且 监 测 精度 达到 毫米 级 别 。 其 原理 如 
Rue) : 

设 研究 区 在 时 间 段 %~ 刀 内 有 V+l 幅 影 像 ,通过 
基线 组 合 得 到 必 景 干涉 图 像 。 


N+1l ye NAID 


Wt, 和 4 时 刻 (b5>a) 得 到 的 图 像 干涉 处 理 后 生 
成 第 i 幅 干 涉 图 , 像 元 (x, y) 处 的 干涉 相位 表示 为 : 
8y) =P, 7P, = pt,%,y) -9 (tix,y) ~ 
Anfd(t,,«.y) = d(t,,x.y VA 
WRP :d(t,.«, y) Fl d(t,, x, y AMS FE SMAI ty 
t, 和 时间 点 对 应 LOS 向 地 表 形 变量 ; p, 和 9g, 为 


(3) 


对 应 时 间 的 解 缠 相 位 ;4 为 雷达 波长 。 

解 缠 后 的 干涉 图 相位 为 : 

6,(x,y) =A (x,y) (4) 

式 中 : 6, 为 干涉 相位 ; A, 为 系数 矩阵 ;@，, 力 为 像 元 
位 置 。 

本 研究 利用 SBAS-InSAR 技术 对 Sentinel- 1 Ft 
轨 数 据 进 行 处 理 并 获取 形变 结果 ,数据 处 理 过 程 如 
图 3 所 示 。 
2.3 精度 评价 

以 20 m 为 距离 对 1990 一 2020 年 冰川 边界 做 缓 
冲 区 处 理 , 用 缓冲 区 面积 除 以 冰川 总 面积 得 到 误差 
率 ( 表 3)。 最 终 得 到 的 误差 率 较 小 ,可 以 满足 精度 
BRO), 


图 2 1990 一 2020 年 部 分 区 域 冰 川 变化 
Fig. 2 Part of glacier variations from 1990 to 2020 


数据 


Fr 一 ”| 相位 解 缠 < | 最 小 费用 流 


ý 
高 程 误差 线性 形变 相位 


y 


残余 相位 解 缠 


生成 连接 图 
/ 参考 DEM 人 > 差分 干涉 处 理 
提取 相干 点 


非 线性 形变 相位 


SVD 估 计 形 变 


ee 


时 间 序 列 形变 


EE: DEM 为 数字 高 程 模型 ;SVD 为 奇异 值 分 解法 。 
图 3 差分 干涉 测量 短 基 线 时 序 分 析 (SBAS-InSAR ) 技 术 流 程 图 
Fig.3 Workflow diagram of SBAS-InSAR technology processing 
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表 3 1990—2020 年 珠 峰 自然 保护 区 缓冲 区 
面积 及 误差 率 
Tab.3 Buffer area and error rate of Qomolangma 
Nature Reserve from 1990 to 2020 


年 份 总 面积 km 缓冲 区 /km 误差 率 /% 
1990 1306.31 78.75 6.03 
2000 1294.82 80.49 6.22 
2010 1231.31 69.67 5.66 
2020 1059.15 72.24 6.82 


3 结果 与 分 析 


3.1 冰川 总 面积 分 布 与 变化 特征 

1990—2020 年 珠 峰 自然 保护 区 冰川 面积 整体 
呈 退 缩 状 态 ,由 1306.31 km 减少 到 1059.15 km’, 共 
退缩 247.16 km’, 面积 变化 率 为 -18.92%。 不 同时 期 
的 冰川 退缩 程度 有 所 差异 ( 表 4)。 其 中 ,1990 一 
2000 年 冰川 面积 退缩 了 11.48 km’, 年 均 变 化 率 为 
-0.09%a';2000 一 2010 年 冰川 面积 退缩 63.51 km ,年 
均 谈 化 率 达 -0.49%…a';2010 一 2020 年 冰川 面积 退 
缩 了 172.17 km?, 相 较 于 前 2 个 时 段 ,该 时 段 冰川 退 
缩 速率 显著 升 高 。 综 上 , 近 30 a 来 , 珠 峰 自然 保护 区 
冰川 旦 加 速 退缩 趋势 , 且 近 10 a 退缩 趋势 愈 发 显著 ， 
与 喜马拉雅 山地 区 冰川 退缩 趋势 一 致 交 。 


表 4 1990 一 2020 年 珠 峰 自然 保护 区 冰川 
面积 分 布 和 变化 
Tab.4 Distribution and change of glacier area in 
Qomolangma Nature Reserve from 1990 to 2020 


qg SERO E 面积 生化 率 ”年 均 变 化 率 
/km /km 1% 1%ora' 
1990 1306.31 7 E 一 
2000 1294.82 -11.48 —0.88 —0.09 
2010 1231.31 -63.51 -4.90 —0.49 
2020 1059.15 -172.17 -13.98 -1.40 
总 计 Š -247.16 -18.92 —0.63 


注 :“-" 表 示 无 数据 。 


3.2 不 同 海拔 的 冰川 分 布 与 变化 特征 

结合 DEM 数据 ,以 200 m 为 间隔 ,统计 1990 一 
2020 年 珠 峰 自然 保护 区 冰川 面积 在 随 海拔 的 分 布 
情况 (图 4a)。 随 着 海拔 升 高 ,冰川 面积 分 布 均 呈现 
先 增 后 减 , 且 各 时 期 分 布 面积 均 在 5600~5800 m 处 
达到 最 大 。 以 1990 年 为 例 ,海拔 3600~4800 m 和 
7000~8200 m 处 冰川 | 分布 面 积 较 少 ,分 别 为 14.35 km? 


和 12.77 km?;5600~6000 m 范 围 内 面积 分 布 较 多 , 达 
445.51 km , 占 冰 川 总 面积 的 34.10%。 海 拔 4800 m 
以 下 冰川 面积 分 布 较 少 ,可 能 是 由 于 温度 影响 且 水 
汽 补给 较 少 ,导致 低 海拔 呈 干 冷气 候 , 不 利于 冰川 
的 发 育 与 积累 ;5200~6400 m 处 有 良好 的 空间 位 置 
和 充沛 的 降水 ,气温 随 海 拔 升 高 而 降低 ,为 冰川 的 
发 育 提供 适宜 条 件 。 冰 川 界 认为 一 座 山 峰 可 能 存 
在 一 个 或 者 更 多 降水 带 ”。 高 海拔 山地 冰川 的 降 
水 量 随 高 度 变化 有 明显 变化 规律 , 随 着 海拔 升 高 ， 
山坡 的 湿润 空气 温度 下 降 ,饱和 度 增加 , 当 到 达 一 
定 高 度 以 降水 形式 下 落 ,形成 第 一 个 降雨 带 ,高 海 
拔 地 区 的 冰川 和 积 雪 具有 较 低 温度 和 充足 的 水 汽 ， 
可 降低 过 境 云 团 温度 ,从 而 形成 降水 , 即 第 二 大 降 
水 带 。 谢 自 楚 等 ”认为 受 局 部 环流 影响 ,珠穆朗玛 
峰 地 区 存在 较 明显 的 “第 二 大 降水 带 ”。 日 本 学 者 
Yasunari 等 3 对 喜马拉雅 山 研 究 发 现 第 二 大 降水 带 
分 布 在 海拔 5000 m 以 上 。 因 此 ,5200~6400 m 范 围 
的 冰川 积累 可 能 是 依赖 于 第 二 大 降水 带 的 补给 。 

1990—2020 年 珠 峰 地 区 冰川 面积 总 退缩 量变 
化 表现 为 先 增 后 减 (图 4b) ,在 5400~$600 m 处 达到 
最 大 ,为 50.06 km。 海 拔 8000m 以 上 几乎 无 退缩 。 
相 比 于 退缩 量 ,研究 区 内 冰川 总 退缩 率 则 呈 波 动 下 
降 态势 ,在 5000~5200 m 处 出 现 小 高 峰 , 达 30.93% 。 
可 以 发 现 : 低 海 拔 区 冰川 退缩 量 较 少 ,但 退缩 率 较 
高 ; 而 高 海拔 处 冰川 退缩 量 和 退缩 率 较 低 。 可 能 是 
因为 低 海拔 处 温度 较 高 ,不 利于 冰川 积累 ,导致 冰 
川 虽然 退缩 量 小 ,但 在 较 少 冰川 分 布 面积 影响 下 ， 
形成 高 退缩 率 。 而 高 海拔 区 温度 相对 较 低 ,冰川 不 
易 消融 ,但 积累 空间 较 小 ,因此 退缩 率 和 退缩 量 均 
不 明显 。 
3.3 不 同 坡度 的 冰川 分 布 和 变化 特征 

以 5° 为 间隔 ,将 冰川 坡度 分 为 12 个 等 级 。4 个 
时 期 冰川 面积 分 布 特征 相似 , 随 坡 度 增加 ,冰川 分 
布 面积 均 表 现 为 先 增加 后 减 小 趋势 (图 5a)。 以 
1990 年 为 例 , 冰 川 在 10°~15° 面 积分 布 最 大 , 达 
196.21 km’, 占 对 应 时 段 冰 川 总 面积 的 14.94% ,55?~ 
60° 分 布 最 少 , 仅 为 21.11 km ,面积 占 比 为 1.62% 。 
可 见 , 珠 峰 自然 保护 区 冰川 主要 分 布 在 较 平缓 地 
带 , 极 平缓 和 陡峭 区 域 分 布 较 少 。 近 30 a 来 研究 区 
冰川 总 退缩 量 和 总 退缩 率 均 呈 现 单 峰 形态 (图 5b)， 
在 10°~15° 坡 度 范 围 内 退缩 面积 最 大 , 达 37.95 km’, 
55°~60° 范 围 内 退缩 量 较 小 , 仅 为 2.53 km*。 而 退缩 
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Fig. 4 Distribution and variation of glacier area at different altitudes from 1990 to 2020 
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图 5 1990 一 2020 年 珠 峰 自然 保护 


(b) 不 同 坡度 面积 变化 特征 
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区 不 同 坡度 冰川 分 布 特征 和 变化 特征 


Fig. 5 Distribution and variation characteristics of glacier area at different slopes 


in Qomolangma Nature Reserve from 1990 to 2020 


率 在 20>~25° 处 最 大 ,55°~60? 范 围 内 最 小 。 
3.4 冰川 形变 特征 

基于 SBAS-InSAR 技术 获取 2020 年 1 月 1 日 至 
12 月 26 日 珠 峰 自然 保护 区 冰川 表面 形变 信息 , 结 
表明 2020 年 珠 峰 自然 保护 区 冰川 平均 形变 速率 介 
于 -129.069~140.252 mma Z EJ (16). 

大 部 分 冰川 形变 呈现 下 降 , 仅 少 部 分 区 域 表 
现 为 抬升 。 选 取 a、.b、c 和 d 4 个 区 域 ,分 析 可 知 :a 
区 域 形 变速 率 达 -129.069~113.944 mm- a`, b K EÈ 
形变 速率 为 -32.811~36.114 mm: a!,c 区 域 形 变 


速率 为 -37.397~36.357 mm-a”, d KREE ERX H 
-21.033~32.013 mm*a'。 可 见 , 大 部 分 冰 舌 区 域 呈 
减 薄 趋 势 ,而 高 海 poe as. 

由 于 缺乏 实地 测量 数据 ,本 研究 利用 非 冰 川 区 
形变 结果 对 冰川 形变 结果 进行 间接 验证 ,在 非 冰 川 
区 域 选 取 4 个 裸 地 区 (图 7) ,统计 形变 量 ( 表 5)。 计 
算得 到 均 方 根 误差 分 别 为 0.16 mm,0.16 mm,0.17 mm 
和 0.17 mm, 证 明 形变 结果 具有 较 高 的 可 靠 性 。 

3.4.1 不 同 海拔 冰川 表面 形变 特征 基于 ArcGIS 软 
件 ,利用 分 区 统计 工具 获取 不 同 海拔 的 冰川 形变 量 
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注 : 工 工 . 开 为 非 冰 川 区 。 
图 7 裸 地 区 域 选 取 


Fig. 7 Selection of bare ground areas 


区 


(图 8)。 结 果 表 明 :不 同 高 程 带 冰川 晨 现 出 不 同形 
变 趋 势 , 人 研究 区 内 绝 大 部 分 冰川 呈现 减 薄 态势 。 海 
IR 5800 m 以 下 冰川 表现 不 同 程度 沉降 趋势 ,其 中 
4200~4400 m 处 沉降 最 为 严重 ,变化 量 达 -17.246 mm. 
海拔 5800~6400 证 出 现 轻微 抬升 ,并 在 6200~6400 m 
处 抬升 较 大 ,形变 量 为 0.347 mm。 上 文中 提 及 珠 峰 
地 区 在 5200~6400 m 处 可 能 存在 “第 二 大 降水 带 ”。 

根据 冰川 形变 结果 分 析 , 冰 川 在 5800~6400 mm 呈现 


表 5 裸 地 区 域 形变 统计 


Tab. 5 Deformation statistics of bare land area 


区 域 选 取 平均 值 /mm 标准 差 /mm ” 均 方 根 误差 /mm 
I 1.58 0.77 0.16 
I =1.25 0.77 0.16 
亚 3.42 0.78 0.17 
TV -0.50 0.78 0.17 
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图 8 2020 年 珠 峰 自然 保护 区 不 同 海拔 冰川 形变 特征 
Fig. 8 Characteristics of glacier deformation at different 


elevations in Qomolangma Nature Reserve in 2020 


抬升 ,结合 前 文 对 冰川 面积 的 探讨 ,由 此 推测 , 珠 峰 
自然 保护 区 的 “第 二 大 降水 带 " 大 致 位 于 $800~6400 m 
处 ,与 Yasunari 等 ”在 喜马拉雅 山 发 现 “ 第 二 大 降水 
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带 ” 的 观测 事实 相符 。6400 m 以 上 冰川 沉降 量 随 海 
拨 升 高 波动 增加 ,在 8000~8200 m 处 形 谈 量 达 
-12.448 mm。 低 海拔 区 域 受 气温 较 高 .湿度 较 小 等 
因素 影响 ,冰川 出 现 消 融 ,导致 冰川 减 薄 。 而 通常 
高 海拔 地 区 由 于 空间 充足 ,降水 充沛 ,冰川 更 容易 
积累 ,从 而 表现 抬升 。 但 本 研究 中 高 海拔 区 域 冰川 
沉降 严重 ,推测 可 能 是 研究 区 发 生 过 雪崩 导致 形变 
异常 。 经 查证 ,研究 区 包含 的 希 夏 邦 马 峰 地 区 曾 于 
2020 年 8 月 发 生 过 较 大 雪 前 。 总 体 而 言 ,冰川 在 低 
海拔 地 区 更 容易 出 现 物质 消融 导致 冰川 减 薄 ,而 高 
海拔 多 呈现 抬升 趋势 。 

3.4.2 不 同 坡度 的 冰川 表面 形变 将 征 以 5° 为 间 
隔 , 将 冰川 分 为 13 个 等 级 (图 9), 珠 峰 自 然 保护 区 冰 
川 在 各 个 坡度 呈现 不 同形 变 程度 。 坡 度 0~50° 内 ， 
研究 区 冰川 旦 现 沉降 态势 ,形变 量 随 坡度 增 大 呈 
“ 减 - 增 - 减 - 增 " 趋 势 , 其 中 40°~45° 范 围 内 ,沉降 最 
为 严重 ,形变 量 达 -1.401 mm。 当 坡度 大 于 50°? 后 变 
化 趋势 减 小 ,在 坡度 5$" 左 右 开始 出 现 抬升 趋势 ,并 
于 大 于 60?" 处 达到 形变 量 最 大 值 , 为 1.545 mm. H 
现 上 述 情况 的 可 能 原因 是 :陡坡 一 般 位 于 高 海拔 
区 ,温度 较 低 ,冰川 消融 速率 较 慢 ,同时 由 于 积 雪 的 
作用 ,导致 冰川 积累 ,抬升 趋势 显著 。 
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图 9 2020 年 珠 峰 自然 保护 区 不 同 坡度 冰川 形变 特征 


Fig. 9 Characteristics of glacier deformation at different 


slopes in Qomolangma Nature Reserve in 2020 


4 讨论 


冰川 面积 变化 对 气候 的 响应 具有 一 定 沾 后 性 ， 
但 气候 变化 规律 可 以 透 过 冰川 物质 平衡 直观 呈现 


出 来 。 近 几 十 年 来 青藏 高 原 气 温 持 续 上 升 , 且 增 温 
速率 大 于 北半球 中 纬度 区 域 近 。 但 近 10 a 青藏 高 原 
东部 增 温 趋势 开始 减缓 ,而 中 西部 地 表 增 温 较 为 强 
烈 。 任 机 文 等 “研究 表明 喜马拉雅 山中 段 气 温 
不 断 波动 上 升 且 在 夏季 升温 更 为 明显 。 过 去 几 十 
年 珠 峰 地 区 的 增 温 速 率 也 显著 大 于 中 国 和 世界 的 
平均 水 平 ”。 青 藏 高 原 降水 变化 受 西风 、 南 亚 风 及 
高 原 季 风 影 响 ,喜马拉雅 山中 段 降 水 主要 源 于 印度 
季风 ,但 研究 区 绝 大 部 分 冰川 位 于 喜马拉雅 山北 
坡 ,由 于 高 大 的 山体 阻挡 了 季风 降水 ,导致 北 坡 降 
水 减少 。 此 外 , 冰 芯 积累 量 可 以 直接 反映 降水 量变 
化 。 侯 书 贵 等 基于 远东 绒布 冰川 冰 世 记录 指出 
珠 峰 地 区 冰川 积累 量 不 断 减少 ,与 冰川 退缩 趋势 一 
致 。 段 克勤 等 ”根据 达 索 普 冰 忌 积累 量 研究 推测 ， 
喜马拉雅 山中 部 地 区 的 季风 降水 将 随 气温 上 升 而 
减少 。 本 研究 中 , 珠 峰 自 然 保护 区 冰川 面积 在 近 30 a 
来 不 断 退缩 , 且 2020 年 冰川 整体 出 现 减 薄 , 物 质 平 
衡 呈 亏损 态势 。 其 呈现 的 物质 平衡 亏损 状态 与 前 
人 研究 一 致 好 ,上 且 从 时 间 尺 度 看 , 珠 峰 地 区 冰川 消 
融 在 进一步 加 强 。 研 究 区 冰川 变化 趋势 与 近年 气 
温 升 高 和 降水 减少 具有 一 定 相关 性 ,因此 ,气温 增 
加 、 降 水 减少 可 能 是 导致 冰川 退缩 的 主要 因素 。 

除 气温 和 降水 影响 外 ,冰川 所 处 地 理 位 置 、 地 
形 差 异 等 因素 ,会 导致 冰川 接收 到 的 太阳 辐射 .水 
汽 来 源 和 季风 强度 不 同 ,进而 影响 冰川 变化 。 由 于 
人 研究 区 降水 受 印 度 季 风 影 响 ,研究 区 范围 较 大 ,而 
季风 降水 具有 空间 差异 性 ,在 一 定 程 度 上 会 导致 
物质 平衡 差异 。 此 外 ,冰川 退缩 与 冰 湖 扩张 之 间 存 
在 密切 联系 , 珠 峰 自然 保护 区 冰川 加 速 退缩 , 冰 湖 
迅速 扩张 中 。 温 度 升 高 导致 冰川 出 现 一 定 程度 消 
融 , 冰 川 末 端 表现 尤为 显著 ,而 末端 冰川 消融 会 导 
致 积 雪 融 水 增加 ,促进 末端 冰 湖 扩张 。 同 时 ,温度 
上 升 导致 冰 湖 蒸发 量 增 加 , 当 莹 发 量 远 小 于 融 水 的 
补给 量 , 冰 湖 继续 扩张 。 可 见 ,冰川 消融 对 冰 湖 扩 
张 有 显著 影响 ,同时 冰川 融 水 也 是 冰 湖 补给 的 重要 
FUR; 而 冰 湖 不 断 扩 张 对 冰川 末端 的 消融 产生 正 
向 作用 ,导致 冰川 末端 加 速 融 化 。 

珠 峰 自然 保护 区 内 分 布 大 量 表 磺 覆盖 型 冰川 ， 
相 比 裸 冰 区 域 , 表 矿 覆盖 型 冰川 表面 反射 率 较 低 、 
吸 热 快 ,可 向 表 厂 下 部 冰川 和 周围 冰川 传递 较 多 热 
量 , 促 进 冰川 消融 。 厚 度 较 薄 的 表 厂 覆盖 对 冰川 消 
融 有 一 定 促进 作用 ,大 于 一 定 厚度 时 有 隔 热 作 用 ” 。 
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表面 , 当 冰 川 表面 污染 , 表 奔 覆盖 面积 进一步 扩张 ， 
冰川 反照 率 继续 下 降 , 吸 热能 力 增 强 , 加 速 冰川 
消融 。 

综 上 所 述 ,气温 升 高 .降水 减少 可 能 是 导致 珠 
峰 自 然 保 护 区 大 部 分 冰川 物质 亏损 的 主要 因素 。 
同时 ,气候 空间 差异 .地形 等 因素 也 可 能 是 导致 物 
质 平衡 差异 的 重要 因素 。 


5 结论 


本 文采 用 多 源 琐 感 数据 对 珠 峰 自然 保护 区 冰 
川 分 布 和 变化 特征 进行 研究 ,并 分 析 冰 川 对 气候 变 
化 的 响应 关系 ,结论 如 下 : 

(1) 1990 一 2020 年 珠 峰 自然 保护 区 内 冰川 面积 
持续 退缩 ,上 且 近 10 a 来 退缩 趋势 愈 发 显著 , 近 30 a 来 
研究 区 内 冰川 总 面积 共 退 缩 247.16 km ,年 均 变化 
率 为 -0.63%a'。 

(2) 珠 峰 自然 保护 区 冰川 主要 分 布 于 海拔 5400~ 
6200 m 和 坡度 10°?~15° 处 ,分 别 占 对 应 时 期 冰川 总 
面积 的 34.10% 和 14.94%, 其 中 5400~5600 m 和 10?~ 
15° 退 缩 最 为 显著 ,分 别 退 缩 50.06 km 和 37.95 km’, 

(3) 珠 峰 自然 保护 区 2020 年 平均 形变 速率 在 
-129.069~140.252 mm a :之 间 , 其 中 海拔 4200~ 
4400 m Fil 40°~45° Mh ,沉降 最 为 严重 ,形变 量 分别 
为 -17.246 mm 和 -1.401 mm. 

(4) 气温 上 升 、 降 水 减少 可 能 是 导致 珠 峰 自然 
保护 区 大 部 分 冰川 物质 亏损 的 主要 因素 。 同 时 , 气 
候 空 间 差异 、 地 形 等 因素 也 可 能 导致 物质 平衡 差异 
的 重要 因素 。 
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Glacier monitoring in Qomolangma Nature Reserve based 


on multi-source remote sensing data 


JI Qin, ZHANG Cuilan, DING Yuekai, CAO Xianggin, LIANG Wenli 
(Chongqing Key Laboratory of GIS Application, Chongqing Normal University, Chongqing 401331, China) 


Abstract: Glacier mass balance is a crucial indicator of climate change and is of great significance for assessing 
the regional ecological environment, thereby preventing and controlling glacier disasters. Based on Landsat series 
images, the ratio threshold method and visual interpretation method are applied to extract the glacier boundaries 
of Qomolangma Nature Reserve from 1990 to 2020. Moreover, the distribution and change characteristics of the 
glacier area in the past 30 years are investigated while the regional glacier deformation characteristics are moni- 
tored based on SBAS-InSAR technology to invert the changes in the glacier mass balance. The following results 
were observed. (1) From 1990 to 2020, the glacier area in the Qomolangma Nature Reserve continuously retreat- 
ed, with this trend becoming much more prevalent in the last 10 years. Moreover, the total glacier area shrank by 
247.16 km’ with a change rate of —18.92%. (2) The glaciers in the Qomolangma Nature Reserve were mostly situ- 
ated at an altitude of 5400-6200 m and a slope of 10°-15°, and the highest ice loss occurred at an altitude of 
5400-5600 m and a slope of 10°-15°. (3) In 2020, the average glacier deformation rate was between —129.069 mm .a ' 
and 140.252 mm -a '. The subsidence and surface deformation of glaciers are most severe at altitudes of 4200- 
4400 m and a slope of 40°—45°. (4) Rising temperatures and decreasing precipitation are believed to be the main 
causes of most glacier material losses in the Qomolangma Nature Reserve. Meanwhile, spatial climate and topo- 
graphic differences may affect mass balance changes. 

Key words: multi-source remote sensing; SBAS-InSAR; glaciers change; mass balance; Qomolangma Nature 
Reserve 


